A high-velocity transient outflow in 77 Carinae 
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We analyze velocity profiles of the X-ray spectral lines emitted by the 77 Cari- 
nae (77 Car) stellar binary at four epochs, just before the X-ray minimum (associ- 
ated with periastron) and more than two years before the minimum (~ apastron). 
The profiles are nicely resolved by the HETGS spectrometer on board Chandra. 
Far from periastron, we find symmetrical lines that are more or less centered 
at zero velocity. Closer to periastron, the lines broaden, shift towards the blue, 
and become visibly asymmetric. While the quiescent X-ray emission and slight 
(< 200 km s _1 ) centroid shifts can be ascribed to the ordinary continuous binary 
wind interaction and to the orbital velocity of the secondary star, the observed 
high-velocity emission up to ~2,000 km s _1 and the abrupt flares during which 
it occurs can not. This leads us to interpret the high- velocity flaring emission as 
due to a fast collimated outflow of ionized gas. 



Subject headings: techniques: spectroscopic 
- stars: individual: 77 Carinae 



X-rays: stars — stars: mass loss 



1. INTRODUCTION 

Despite being studied for many yea rs and in all waveba nds, the 77 Carinae system (7/ Car), 
best known for its spect acular nebula (IMorse et al.lll998f ) that can be tracked back to the 
Great Eruption of 1840 ( Davidson &: Humphreys 19971 ) is still much of a mystery. It is now 
believed to be a stellar binary system with a period of 5.54 years (jDaminelil Il996l ). The 
primary star is one of the brightest and most massive (~ 120 M ) stars in our galaxy. 
Although the X-ra ys detected from 77 Ca r can be largely explained by the collision of the 
two stellar winds (j Corcoran et al.l l2001al : iPittard fc Corcoran! 120021 ) , the ~ 70 day X-ray 
minimum that is periodically observed around periastron is still not well-understood. 



Currently, the mass loss rate of the primary exceeds 1 4 M^ vr 1 rtSmith et al.l 120031 ) 



and that of the secondary is an order of magnitude lower ( Pittard fc Corcoran 20021 ). The 
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wind of the secondary is much faster (few 1,000 km s 1 . |Pittard fc Corcoran! |2002|) than that 
of the primary (~ 600 km s _1 and depending on latitude. ISmith et al.ll2003l ). The X-ray flux 



between 2 and 10 keV during the past 12 years is normall y at a level of 5x10 erg s 



-1 —2 

cm , 



as expected from the coll i sion of the two ste llar winds (lUsovl Il992l ; ICorcoran et al.l l2001al ; 
Pittard fc Corcoran! |2002| ; lAkashi et al.ll2006l ). The X-ray light curve is roughly constant 
throughout most of the orbit, but rises gra dually by a factor o f 3-4 upon approach to 
periastron at which time short flares appear (j Corcoran et al.lll997l ). before it drops sharply 
and stays in an X-ray low state (~ 10% of normal bright ness) for approximately 70 days. The 
full X-ray light curve can be found in ICorcoran! (120051 ). Although absorption by the dense 
primary wind plays a certain role in attenuating the X-ray emission from rj Car, both spectral 
and temporal arguments have been put f orward to reject absorption or an eclipse as the 



prim e reason for this persisting low state (jHamaguchi et al.l 12004 ; ISokerll2005l ; lAkashi et al. 



20061 ). which is unique to rj Car. It has been proposed that if the binary separation during 



( Soker 


2005; 


Akashi et al. 


2006) 



of the broad-ban d X-ra y behavior around the 2003 periastron passage was conducted by 



Hamaguchi et al.l (120071 ). who claim the unchanged temperature during the X-ray minimum 
argues against accretion. However, i f the X-ray emissio n during the minimum is due to the 
secondary's residual fast polar wind (j Akashi et al. 2006 ). it would indeed result in unchanged 
X-ray temperatures. 

In this work, we wish to add to the temporal aspects of the rj Car tale a high-resolution 
X-ray spectroscopic dimension. This is achieved by using five deep archival exposures of 
r) Car with the High Energy Grating Spectrometer (HETGS) on board the Chandra X- ray 
Observatory (CXO), only the first of which has been published ( Corcoran et al. 2001al lbh . 
Theoretical variations of spect ral line profiles in binary colliding wind systems have been 
modeled by lHenley et al.l (120031 ) . The varying line profiles presented here for r] Car, however, 
do not follow these models, at least not in a straightforward way, as we explain in §3. 



2. OBSERVATIONS AND RESULTS 

The log of the Chandra observations is given in Table 1. Two observations (3745 and 
3748) coincide with the short intense flares just before periastron. The exact orbital param- 
eters of r] Car are not well determined, but we assume throughout this paper a period of 
5.54 years (2024 days) and periastron to be at 2003-06-29. The X-ray low state, thus, occurs 
during = - 0.035. The data were retrieved from the CXO archive and processed using 
standard software (CIAO version 3.2.2.). Chandra 1 ?, HEG and MEG (high and medium en- 
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ergy gratings) spectrometers operate simultaneously in the 3.5 - 7.5 A band of interest here. 
Aiming at the highest possible kinematic resolution, the higher spectral resolving power of 
the HEG in this band (factor of ~ 2) is favored over the slightly higher effective area of the 
MEG (factor of ~ 2), especially since the co-added plus and minus 1st order HEG data are 
already of sufficiently high signal-to-noise ratio (S/N). 

To reveal the dynamics of the X-ray gas, we transferred each spectral line profile to 
line-of-sight velocity space (vu in km s _1 ) using the non-relativistic Doppler shift vn = 
c(A — A )/A , where A and A represent the observed and rest-frame wavelengths of the 
line, respectively, and c is the speed of light. Negative and positive velocities represent ap- 
proaching and receding gas, respectively. Fig. 1 shows the broad and variable X-ray line 
profiles of r) Car demonstrated on the Si +13 Lyct doublet at 6.18 A. As expected, the MEG 
profile represents a smoothed version of the higher-resolution HEG profile. The observed 
line is much broader than the instrumental line spread function (also shown). Fig. 1 also 
demonstrates how the line profile varies dramatically from phase = -0.470 (Obs. 632) to 
= -0.028 (Obs. 3745). The early-phase (0 = -0.470) line is centered around zero velocity 
with no significant emission beyond 700 km s _1 . The line closer to periastron (0 = -0.028) 
features high velocity gas up to ~2000 km s _1 . 

In each observation, the kinematic profiles of the bright spectral lines are fairly similar. 
This is demonstrated in Fig. 2, where the four best S/N lines of Obs. 3745 (0 = -0.028) are 
shown. In order to elucidate how the line profiles vary as a function of phase, we co-added 
the profiles of nine bright lines in each observation, namely the Lya (ls-2p), Hea resonance 
(ls 2 -ls2p) and Hea forbidden (ls 2 -ls2s) lines of the Si, S, and Ar ions. The similarity 



of the profiles of these lines in a given observation (see also iHenleyl 120051 . Figs. 4.9-4.12) 
justifies the co-adding procedure, which yields an average profile with improved S/N. The 
sufficiently large separation of the resonance to forbidden lines (3400 - 4100 km s _1 ) and the 
weak intensity of the intercombination lines in the He-like triplets ensure that the average 
profiles are not confused by blends. Individual Fe lines below 2 A could not be included in 
this analysis, despite their high intensity, since the HETGS resolving power decreases strongly 
at these low wavelengths to the extent that the He-like Fe line complex is unresolved. The 
blue wing of the Hea Fe resonance line is consistent with, though not resolved as well as, the 
other line profiles presented here. The resulting mean profiles have been normalized (except 
for Obs. 3747) to facilitate the comparison and are presented in Fig. 3. 

A clear trend can be seen in Fig. 3. Far from periastron, the line profile is relatively 
symmetric and narrow. Closer to periastron, bright components of gas moving towards us 
at velocities as high as -2,000 km s _1 start to develop. As the system further approaches 
periastron, the outflow dominates the line profile to the point where the bulk of the emission 
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is clearly blue-shifted. After periastron, the emission lines are consid erably weaker. Lin e 
profiles from a binary wind system are not expected to be Gaussian (IHenley et al.l 120031 ) . 
Nevertheless, in order to get an idea of the line shifts and widths, we tried to fit Gaussians to 
the mean profiles in the -2,500 km s _1 to +1,300 km s _1 range (Fig. 3). The results are pre- 
sented in Table 2. For the two late-phase observations (0 = -0.028, -0.006) single Gaussians 
provide poor fits that can be significantly improved by adding a second component. This 
does not imply there are necessarily exactly two velocity components, but merely that the 
velocity distribution at these phases is more complex than Gaussian. The double-Gaussian 
best fits are plotted in Fig. 3. A second isolated peak at (ft = -0.028 is not as conspicuous 
as at = -0.006. Hence the broad high-velocity component and the larger errors at <ft = 
-0.028 (Table 2). 



3. INTERPRETATION AND DISCUSSION 

The colliding wind binary model, which is generally accepted for explaining the X-ray 
emission of rj Car, has several clear predictions regarding the spectral line profiles expected 
at different orbital phases. The Doppler shifts are the result of shocked gas flowing along the 
contact discontinuity (CD) surface and away from the stagnation point (SP). If the shock 
opening half- angle is between 45° — 90°, the broadest line profiles are expected when the 
system is observed perpendicular to the line connecting the two stars, as the hottest gas 
formed around the SP flows directly towards and away from the observer. In this case, the 
lines are centered at zero velocity and are intrinsically symmetrical, but could be skewed 
by absorption of the far (red) emission. On the other hand, when the system is observed 
from behind one of the stars, shifted centroids are expected. If the system is observed from 
behind the weaker (stronger) wind, downstream shocked gas flows along the CD surface in 
the general direction of (opposite) the observer, producing blue- (red-) shifted centroids. The 
maximum velocities in this case would be less than those in the first (orthogonal ) orientation. 



For a complete and detailed calculation of line profiles from colliding winds see lHenley et al. 



( 120031 ). In the following, however, we argue that the current line profiles observed at the 
different phases of the r] Car orbit are in contradiction, even qualitatively, with the simple 
colliding wind scenario. 



If the major axis of the 77 Car binary is perpendicular to our line of sight (ISmith et al. 



20041 ). the broadest profiles are expected at apastron (see above), but that is when the 



observed profiles are narrowest (<p = -0.470). Also, with this geometry, the secon dary would 



have to pass in front of the primary before periastron (c.f., ISmith et al.l 120041 ) to explain 



the blue-shifted centroids at = -0.028, -0.006, but that would produce symmetric (blue- 
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shifted) profiles as absorption through the secondary wind is weak, while asymmetric profiles 
are observed. In any case, the lower velocities observed near apastron appear to rule out 
this orientation if the colliding wind geometry is to produce the observed X-ray line profiles. 

Alternatively, if the projection of the line of sight on the orbital plane is along the majo r 
axis and the viewing angle at apastron is from behind the secondary (ICorcoran et al.ll2001al ). 
the observed blue-shifts at <fi = —0.028, —0.006 could be ascribed to gas flowing from the SP 
towards the observer and the unobserved red wing of the line to absorption by the primary 
wind. However, according to this scenario, at apastron (0 = —0.470) one would expect 
blue-shifted centroids of at least a few 100 km s -1 due to downstream gas (see above), which 
are not observ ed. For an observer a bove the binary plane (i < 90°), as seems to be the 
case in rj Car (ICorcoran et al.ll2001at ). t he problem becomes m ore severe (higher blue-shift 
is expected, but not observed). Indeed, iKashi fc Sokerl (120071 ) showed that the X-ray line 
shifts with orbital phase can not be due to the change in orientati on of the collid ing winds 
region when the inclination of the orbital plane is considered (c.f., ICorcoran! 120071 ). 



Apparently, the naive wind-collision scenario alone is insufficien t to exp l ain th e observed 
Doppler shifts of the X-ray line profiles of r\ Car at all binary phases. iHenleyl (120051 ) compared 
the same data with line profiles calculated from his hydrody namical simul a tions and reached 
a similar conclusion. A recent attempt has been made by iHenley et al.l (120071 ) to propose 
yet another orientation. However, we note that the asymmetric line profiles are observed 
during strong peaks in the X-ray light curve. These peaks too can not be accounted for by 
the simple colliding wind models. As far as we understand, the high X-ray velocities remain 
unexplained and we ascribe them to a transient outflow that arises just before the onset of 
the X-ray minimum and is coincident with the X-ray flaring. 



Kashi fc Sokerl (120071 ) suggested orbital motion to be important before periastron. This 



could account for small X-ray line shifts, but not for the highest velocities observed. Uti- 
lizing the two-Gaussi an fits in Table 2, however, we can confront the orbital solution of 
Kashi fc Sokerl (120071 ) with the low- velocity fitted component (column 2 in Table 2). We 



assume an eccentricity of e = 0.9, stellar masses of Mi = 12OM , M 2 = 3OM , and an ob- 
server 37° above the orbital plane (i = 53°), whose projection on t he orbital plane i s exac tly 
behind the secondary at periastron. This ge ometry was invoked by IKashi &: Sokerl (120071) to 
explain the shifts of the He I lines (but c.f., ICorcoran et al.ll2001al ; IHenley et al.l 120071 ). The 
calculation shows that the orbital velocity of the X-ray gas is ~ 80% th at of the seco ndary. 
We added -8 km s _1 to the model for the systemic velocity of rj Car (jSmithl 12004 ) . The 
resulting theoretical velocities of the X-ray gas as a function of phase are seen in Fig. 4 to 
formally agree with the observed Doppler shifts of the low- velocity component, except at 
cf) = -0.006, where the model velocity is slightly smaller. Given the many uncertainties of 
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the binary parameters, this discrepancy may not be significant. 



The high- velocity X-ray emission (column 4 in Table 2), on the other hand, is not 
explained by the orbital motion and could be due to a transient collimated flow ejected from 
the immediate vicinity of the binary system. The peaks in the X-ray light curve hint that 
the outflow may be in the form of blobs. The unresolved hard X-ray Chandra images restrict 
the outflow to within ~ 2xl0 16 cm of the center (0.5", assuming the distance of 2.3 kpc). 
The appearance of the same charge states during all phases implies that the temperatures of 
the X-ray gas and line-emitting outflow remain in the range of kT = 2-5 keV. Hence, the 
outflow likely consists of gas shocked by the collision of the winds (as observed throughout 
the orbit) and is boosted near periastron. The widths of the major peaks seen in Figs. 2 and 
3 suggest that the outflow is only moderately collimated to within ~ 30°. The compression 
of shocked gas provides a natural explanation for the short intermittent flares in the X-ray 
light curve. A collimated outflow from the secondary has also been sugge sted to explain 
the enhanced He II A4686 emission before periastron (ISoker fc Beharll2006l ). Interestingly, 
both proper motion and Doppler shifts h ave been measured for high- velocity visib le-light 
knots much further out from the center ( IWalborn et al.l Il978l : iMeaburn et al.lll993l ). The 
optical knots were eje cted more than a hundred years ago along the minor axis of the nebula 
(IMeaburn et al.lll993l ). We would expect the present outflow to be directed along the major 
axis of the nebula although we have no pertinent information on its present direction with 
respect to the system's geometry. As the massive rj Car primary is a short-lived star shedding 
considerable mass, it may be a supernova in the making. At the very least, it highlights the 
next periastron passage on 2009, January 12 as a faithfully scheduled experiment for the 
astrophysically common, but poorly understood phenomenon of collimated flows. 
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Table 1. CXO/HETG Observation Log 



CXO archive 


Start Time 


Exposure 


Assumed Orbital 


2-10 keV Flux 


ID a 




(ks) 


Phase (0) b 


(10 -10 erg s -1 cm -2 ) 


632 


2000-11-19 02:46:40 


90.69 


-0.470 


0.50 


3749 


2002-10-16 08:08:49 


93.96 


-0.130 


0.98 


3745 


2003-05-02 11:56:16 


97.29 


-0.028 


2.2 


3748 


2003-06-16 05:35:28 


100.1 


-0.006 


0.97 


3747 


2003-09-26 22:45:53 


72.16 


+0.044 


0.48 



a One additional observation was carried out during the X-ray minimum (ID 3746 at 2003- 
07-20 01:46:22) for which no useful spectrum could be obtained. 

b Phase is approximated based on the assumption of a 5.54 year (2024 day) orbit, where 
phase zero is set at 2003-06-29 and the X-ray low state occurs during phase - 0.035. 



Table 2. Kinematic parameters in km s 1 obtained from Gaussian fits to the mean line 

profiles presented in Fig. 3 



Assumed Phase (0) 


Vi 


FWHMi 


V2 


FWHM 2 


X 2 /d.o.f a 


-0.470 


-25 ± 25 


710 ± 60 






0.36 


-0.130 


-50 ± 20 


730 ± 40 






0.43 


-0.028 


-110 ± 30 


730 ± 100 


-830 ± 110 


2030 ± 400 


0.39 


-0.006 


-180 ± 30 


620 ± 60 


-1030 ± 30 


840 ± 90 


0.68 



a Double-Gaussian fits to the profiles at = -0.470 and at = -0.130 yield an un- 
constrained second component. Single-Gaussian models for the profiles at = -0.028 
and = -0.006 yield visibly inferior fits with x 2 / d.o.f of 1.16 and 1.02, respectively. 
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SiLya A.6.18 




Line-of-sight velocity (km s" 1 ) 



Fig. 1. — Velocity profiles of the Si +13 Lya unresolved doublet (6.180 and 6.186 A). Upper 
panel shows consistent HEG (higher spectral resolution) and MEG (smoother) data from 
Obs. 3745 {(f) = -0.028). The profile is clearly broadened up to 2000 km s -1 , much beyond 
the instrument line spread functions (dashed lines). Bottom panel shows the HEG profile of 
Obs. 632 (0 = -0.470, dashed line) and Obs. 3745 (0 = -0.028). The scaled up line from 
the early phase is rather symmetrical and much narrower, showing no line emission beyond 
700 km s -1 . 



- 11 - 



t 1 1 1 1 1 1 1 1 1 rq 




— I 1 1 1 1 1 1 1 1 1 1 — 

-3000 -2000 -1000 1000 2000 



Line-of-sight velocity (km s" 1 ) 



Fig. 2. — Profiles of four bright spectral lines from Obs. 3745 (0 = -0.028) showing asym- 
metric profiles with blueshifts of up to ~ — 2, 000 km s~ l . Horizontal bars represent the 
shift and FWHM Gaussian-fitting results by iHenleyl (120051 ) for these lines, which are not 
inconsistent with the present profiles. 



-| — I — I — I — I — I — I — I — I — I — I — I — I — r 

—■—phase -0.47' 
—•—phase -0.130 
phase -0.028 
—▼—phase -0.006 
phase +0.044 
(not normalized) 



i — i — i — i — i — i — i 



i — i — i — i — I — i — i — i — i — I — i — r - 




c 
3 



X 

E 

1) 

13 



o 
Z 



~ l — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i — i — i — i — i — r~ 

-2500 -2000 -1500 -1000 -500 500 1000 



Line-of-Sight Velocity (km s" ) 



Fig. 3. — Mean, normalized, velocity profiles constructed from nine different spectral lines 
for the five observations of rj Car. Only the profile of Obs. 3747 (0 = +0.044) during which 
r) Car was strongly absorbed is not normalized and plotted at its correct scale with respect 
to Obs. 632 ((f) = -0.470). Typical la errors on each data point are 10 - 15 % (top right 
hand side). The double- Gaussian fits to the (p = -0.028 and (ft = -0.006 profiles are also 
plotted (see Table 2). 
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Phase Phase 



Fig. 4. — Several theoretical variables as a function of phase assuming the orbital parameters 
in the text, throughout the orbit (left column) and close to periastron (right column). Top 
row: Orbital separation (in AU) and relative orbital speed of the two stars (in units of 10 
km s" 1 on the left and 100 km s _1 on the right). The angle 9 is the relative direction of 
the two stars as measured from periastron (scale in degrees on the right). Bottom row: 
Projected velocity (in km s" 1 ) on the line of sight, of secondary star (thick line) and of the 
X-ray gas (thin line). Vertical lines mark the measured Doppler shifts of the slow component 
in Table 2. 



